Nature of the Bogoliubov ground state of a weakly interacting Bose gas 
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As is well-known, in Bogoliubov's theory of an interacting Bose gas the ground state of the 
Hamiltonian H = X^k^o k ls f°und by diagonalizing each of the Hamiltonians corresponding 
to a given momentum mode k independently of the Hamiltonians ffvuk) of the remaining modes. 
We argue that this way of diagonalizing H may not be adequate, since the Hilbert spaces where the 
single-mode Hamiltonians are diagonalized are not disjoint, but have the k = in common. A 
number-conserving generalization of Bogoliubov's method is presented where the total Hamiltonian 
H is diagonalized directly. When this is done, the spectrum of excitations changes from a gapless 
one, as predicted by Bogoliubov's method, to one which has a finite gap in the k — > limit. 
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Introduction. Since its inception in 1947, Bogoliubov's 
approach to interacting Bose systems [l[ has been one of 
the most influential theories in condensed matter physics. 
[34l3j Yet, for all its notoriety and popularity within the 
physics community, a key aspect of this theory, having 
to do with the decoupled way in which the Hamiltonian 
is diagonalized, is still not fully understood. Indeed, and 
as is well-known, in the standard formulation of Bogoli- 
ubov's theory, the Hamiltonian H is written as a decou- 
pled sum of contributions from different momenta of the 
form H = X^k^o-^k, where each Hamiltonian de- 
scribes the interaction of bosons in the condensed k = 
state with bosons in the momentum modes ±k, then each 
of the single-mode Hamiltonians is diagonalized sep- 
arately and the ground state (GS) wavefunction of H is 
written as the product of the GS wavefunctions of the 
i/k's. In this letter, we shall argue that, while this way 
of diagonalizing the total Hamiltonian H may seem to be 
valid from the perspective of the standard, number non- 
conserving Bogoliubov's method, where the k = state is 
removed from the Hilbert space and hence the individual 
Hilbert spaces where the Hamiltonians {H^} are diag- 
onalized are disjoint with one another, from a number- 
conserving perspective this diagonalization method may 
not be adequate since the true Hilbert spaces where the 
Hamiltonians {Hk} should be diagonalized all have the 
k = state in common. We then shall discuss a varia- 
tional, number-conserving generalization of Bogoliubov's 
theory in which the k = state is restored into the 
Hilbert space of the interacting gas, and where, instead 
of diagonalizing the Hamiltonians separately, we di- 
agonalize the total Hamiltonian H as a whole. When 
this is done, the spectrum of excitations of the system 
changes from a gapless one, as predicted by the stan- 
dard, number non-conserving Bogoliubov method, to one 
which exhibits a finite gap in the k — > limit. 

Variational formulation of Bogoliubov's theory. We 
shall start by discussing a variational formulation 0, [l4[ 
of Bogoliubov's theory which, historically, has con- 
stituted the basis of the justification of the num- 
ber non-conserving formulation of this method. As 



is well-known, in Bogoliubov's approach one only re- 
tains in the total Hamiltonian H of the system (i) 
kinetic energy terms of the form X)k^o £ka ik ak ' fa) 
Hartree terms ^ k k , ?;(0)a k a k ,akak' /2V, (Hi) Fock terms 

Ek^o ' J C c ) a k a k a o o/^ describing the exchange interac- 
tion between condensed bosons and depleted ones, and 
(iv) pairing terms of the form J^k^o v (^)( a o a o a k a -k + 
a o a o a ic a -k)/2X / . (In the above expressions, = 
h 2 k 2 /2m is the kinetic energy of a boson of mass m and 
wavevector k, i>(k) is the Fourier transform of the interac- 
tion potential between bosons, and V is the volume of the 
system. On the other hand, a k and are creation and 
annihilation operators, respectively). Considering that 
we will be focusing on systems having a fixed number of 
particles N, it is convenient to take the origin of energies 
at the Gross-Pitaevskii value v(0)N(N - 1)/2V. Then it 
can be shown [3] that the Hamiltonian can be written as 
a sum of independent contributions from different values 
of k of the form H ~ Sk^o wnere (throughout this 
paper, for all explicit calculations we shall be using the 
interaction potential v(r) = gS(r), for which v(k) = g): 



+ V S ) /„ + : 
H k = £ka k a k + (2a ( ! ) aoa k ak 
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We now proceed to diagonalize the Hamiltonian by 
considering a hypothetical system where bosons are only 
allowed to be in one of the three single particle states 
with momentum k, or — k. In order to formulate a 
variational approach for the Hamiltonian describing 
such a system, it is sufficient to restrict our attention to 
the Hilbert space fiu spanned by kets \n) of the form: 



|n) = 1^—2?-?,, n, n) 



(al) N ~ 2n (atr («-k)" 



y/(N-2n)\ \jn 



|0>, (2) 



having n bosons with momentum k and momentum — k, 
and N — 2n bosons in the k = state. The general ex- 
pression of the GS wavefunction |^k) of the Hamiltonian 
i?k in this Hilbert space is given by |^>k) = J2n=o C n \n), 
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and it can easily be verified that the expectation value of 
i?k in the state j^k) can be written in the form: 



N/2 

W»k|#k|tfk> = Yl {|C«| 2 n.-k -T-"(-V - 'In) 



v{k) 



V 



v(k) 



+ -^nC*C n - iy /(N -2n + 2)(N - 2n + 1) 

+ ^r(n + l)C* n C n+1 V(JV - 2n - 1)(JV - 2n)}, (3) 

where it is understood that C_i = Cx+(j\r/2) = 0. Assum- 
ing, for simplicity, that the coefficients C„ are real, it fol- 
lows that, for v(k) > 0, the expectation value (ipk\Hk\'>Pk.) 
will be lowered if the coefficients C n have alternating 
positive and negative signs. In this case, the terms on 
the second and third line will be negative, making the 
expectation value lower than what one would obtain if 
products of the form C„C„_i and C„C„+i are positive. 
Bogoliubov's theory corresponds to a variational ansatz 
in which the coefficients C n are assumed to be of the form 
C n = (— Ck)", where the constant Ck > is to be deter- 
mined variationally. The coefficients C n are expected to 
decrease with n, which encodes the fact that the proba- 
bility amplitude of states \n) with a large number n 3> 1 
of bosons having a wavevector ±k will be small. This 
implies that the constant Ck must be less than unity. 

Inserting the variational ansatz C n = (— Ck) n into Eq. 
([3]) , and making use of the approximation \J N(N + I) ~ 
N + | which is valid for iV ^ 1, we can write: 



N/2 



v(k) 



(^kl^kl^k) {( C 2 "K + ^n{N - 2n) 

n=0 



+ 



n(c k ) 2 "- X (iV-2n + |) 



+ ^(n + l) (Ck )^(iV-2,-i)}. 



(4) 



The summations in the above equation can be calculated 
analytically by taking successive derivatives with respect 

to the variable x of the result En=Q x " — 1/(1 ~ x )i van d 
for \x\ < 1 and N — > oo, hence obtaining: 



N/2 

E'< 

71=1 



AT/2 



Cl .ri- ^ ~ (1 - xf ' 

Using these last two results in Eq. (j4)), we obtain: 



(5) 



(Vk|#k|Vk) 



Ck 



1 



Ck£k + v(k)nB(ck - l) , (6) 



where we denote by n B = N/V the density of bosons 
in the system. On the other hand, from the definition 
of \ipk), we easily see that the norm of the wavefunction 

(VkhM ^ given by (V> k |^k) = ETo \ C n? - W ~ 
c k ). Now, if we divide Eq. by this last expression 
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FIG. 1: Coefficients C n of the normalized wavefunction \ip^) 
for N = 100 bosons and k = 0.1. The crosses are the results 
of Bogoliubov's theory, and the circles are the results of the 
exact numerical diagonalization of the Hamiltonian Hk. 



°f (V'klV'k)j we can write for the normalized expectation 
value (-ffk)k = (V'k|-ffk|V ; k)/(V'k|V'k) the following result: 



1 



[sk + u(fc)ns] - v{k)n B 



Ck 



1 



(7) 



Minimization of the above expectation value with respect 

to Ck leads to the quadratic equation c k — 2 ^ t ,(f ) k rtB J c k + 

1 = 0, where we defined £k = £k + v{k)n B . The above 
quadratic equation has two roots, of which only one sat- 
isfies the constraint < Ck < 1 for arbitrary values of £ k. 
This root is given by: [14 1 



Ck 



v(k)n B 



Sk-JSl-v{kf 



(8) 



This result for the constant Ck fully determines the coeff- 
cients C n — {— Ck)™ of the variational ground state l^k) of 
the Hamiltonian H^. The coefficients C n — C n y/l — c k 
of the normalized wavefunction |?/> k ) = \ipv) / y (V'klV'k) 
are plotted (as the crosses) in Fig. [T]for k — 0.1 in dimen- 
sionless units such that k 



Hk/ y/2mv(0)riB- The circles 
in this last figure are the coefficients obtained by direct 
numerical diagonalization of the Hamiltonian iJk- It is 
seen that there is a pretty good agreement between the 
results of our variational method and the exact numerical 
diagonalization for the particular value of k chosen. 

The expectation value of Hk in the ground state 
can readily be found if we use the result (j8]) for Ck in Eq. 
([7]), upon which we obtain: 



<#k)k = ~ 



£k[£k + 2n B v(k)] - e k - n B v(k) . (9) 



The result (9) is exactly what one obtains in the stan- 
dard, number non-conserving Bogoliubov approach for 
the expectation value of a given contribution to the 
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Bogoliubov ground state energy. This agrees with the 
well-known fact that Bogoliubov's theory is a theory in 
which the Hamiltonians are diagonalized indepen- 
dently from one another in essentially disjoint Hilbert 
spaces. Indeed, the quantity in Eq. © is nothing but the 
expectation value (i'k(N)\H k \'tp k (N)) , where |V>ic(iV)) i s 
the normalized ground state of the single momentum 
mode Hamiltonian H k . The reason such a result is ob- 
tained is because in the standard formulation of Bogoli- 
ubov's theory, ao and aj are replaced by the c-number 
y/~N. This implies that the commutators [Hk, -ffk'] van- 
ish identically for k ^ k', which allows the ground state 
wavefunction \^ B {N)) of the total Hamiltonian H to be 
written as a product of the ground state wavefunctions 
for each of the Hamiltonians H k |3J (M here denotes the 
total number of momentum modes kept in the calcula- 
tion, which can eventually be taken to infinity): 

oo oo 

\t> B (N))= C ni C n2 ---C nM \n i; ...;n M ), 
m — o 11 m —0 

with |n i; ...;nM) = n (a ^ (a ^ l )- ( 10 ) 



i=l 



Here, we would like to emphasize that the above expres- 
sion of the wavefunction \^ B {N)) is only consistent with 
the variational constants given in Eq. ([5]) when the k = 
state is removed from the Hilbert space, with ao and aj 
being replaced with y/~N. This means that the ground 
state wavefunction of Eq. (p~0|) above, despite the ap- 
pearance of the contrary, corresponds to a number non- 
conserving approach. A major question that arises is to 
know how the above result will change if we restore the 
k = state to the Hilbert space, and if instead of di- 
agonalizing each of the Hamiltonians H k separately, we 
diagonalize H directly. This will be done next. 

Variational approach for the full Hamiltonian H . We 
now want to generalize the variational treatment of the 
single-mode Hamiltonian H k to treat the full Hamilto- 
nian H = X)k^o of the interacting Bose system. To 
this end, we shall use for ^(iV)) the expression: 



oo oo 



|*(jv)) = £•••£ c ni c n2 ...c n 



m— nj\j— 
M 



x \N - 2^n 8 ;n 1 ; . • . ;n M ), (11) 

where the normalized basis wavefunctions are given by 
(compare with Eq. (JTUJ) of the single-mode theory): 



M 



\N - 2y]m;ni; . . .;n M ) 



K) 



t\JV-2E"i«i 



i=l 



[N -2 ni] 



Note that the GS wavefunction in Eq. (|TT|) is not a simple 
product of GS wavefunctions for the single-mode Hamil- 
tonians H k , and that, even though the expression of these 
single-mode Hamiltonians H k are decoupled and com- 
mute with one another, the presence of all the n^s in 
the number of condensed bosons \N — 2 Ej=i n i] ac ^ s 
like an implicit and rather nontrivial coup ling between 
all these Hamiltonians. One can now show [16| that the 
expectation value (i?k ) of the Hamiltonian H k . in the 
state |^(JV)) is no longer given by Eq. ([7]), but by the 
following expression: 

(H kj ) - -, ^ kj 2 l c kj +v(kj)n B -v(kj)n B \, (13) 

1 C k, ^ L J ' 



where v(kj) is given by: 



«(y-*)(i-^Er 



k^O 



(14) 



If it were not for the term between parenthesis in this 
last equation, the result in Eq. (IT3l) would be perfectly 
identical to the expectation value obtained within the 
single- mode approach, Eq. (J7J. It can be shown [l6} 
that minimization of the trial ground state energy given 
in Eq. (TTBl over the constants Ck leads to a solution of 
the form: 

c k = 1 + C^(P + a) - ^/[l + C^ + cr)] 2 -!. 

(15) 

where the constant a is obtained by solving a non- 
linear self-consistency equation obtained from the min- 
imization procedure, [16| and Cd = 1 — N^/N , with 



E 



k#0 c k 



/(I — c 2 .) the total number of depleted 



bosons. To fix ideas, we shall henceforth consider an 
interacting Bose gas in the dilute limit, and fix the pa- 
rameter n B a 3 to be n B a 3 = 10~ 3 (a here being the scat- 
tering length, which is related to the interaction strength 
g through the relation g — 4nah 2 /to). For this particu- 
lar value of n B a 3 , we find a — 0.39 and Cd — 0.9762. In 
order to show that these values of a and Cd do indeed 
correspond to a lower energy than what one would ob- 
tain by using the coefficients Ck of the single-mode the- 
ory from Eq. ||SJ in Eq. fig]), in Fig. [2] we plot the 
product k 2 ( 1 i'(N)\H k \ , i>(N)), which appears in the eval- 
uation of the ground state energy (the factor k 2 coming 
from the Jacobian in spherical coordinates in three di- 
mensions), as a function of the dimensionless wavevector 
k, for the above two choices of the constants Ck- It can be 
seen from this plot that the solution (fl~5|) with nonzero a 
(solid line) leads to a lower value of the ground state en- 
ergy (fl"3|) than the standard solution of the single-mode 
theory with a = from Eq. (jSJ, (dashed line). Since 
the coefficients in Eqs. © and (fTB"]) are quite different 
from one another, it follows that expectation values of 
observables calculated with the coefficients Ck obtained 
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FIG. 2: Plot of the product k 2 x (*|i/ k |*> vs. fc from Eq. 
(|13[l , using the coefficients Ck (i) from the single-mode theory, 
Eq. ©, (dashed line) and (ii) from Eq. Jl5j with a = 0.39 
and C d = 0.9762 (solid line). 
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FIG. 3: Plot of the depletion of the condensate iVk as a func- 
tion of the dimensionless wavevector k using the coefficients 
Ck (i) from the single-mode theory, Eq. ([8}, (dashed line) and 
(ii) from Eq. ([15]) with a = 0.39 and C d = 0.9762 (solid line). 



by minimizing the expectation value of the full Hamilto- 
nian H will be quite different from the expectation values 
of the same observables calculated using the usual Bo- 
goliubov approximation where, for each value of k, the 
expectation value of the single-mode Hamiltonian is 
minimized. 

As a first example, in Fig. [3J we plot the depletion Ak 
of the condensate A k = (*(7V)|aj i a k |*(^)) = 4/(1-4) 
vs. wavevector k, with the solid line representing the re- 
sult one obtains using the coefficients Ck from Eq. ([T5|). 
and the dashed line representing the results one obtains 
using Eq. ([8]). As it can be seen, the two results are 
qualitatively very different, with diverging like 1/k 2 
as k — > in the standard Bogoliubov theory (which is of 
course unphysical for a system of fixed number of bosons 
A, and leads in one spatial dimension to an infrared 
divergence of the total number of depleted bosons N), 
while, on the contrary, when the improved coefficients Ck 
of Eq. ([15)1 are used, Ak is finite for all values of the 
wavevector k. 

As a second example, we consider the energy to excite 
one boson from the condensate to the single-particle state 
with wavevector k. This is the quantity given by: 



(y\aia k Ha ai\y) 
(y\ala k a a{\y) 



-(V\H\$) (16) 



In the standard formulation of Bogoliubov's theory, 
where ao and aj are replaced by the c-number ^/Nq 
(where N is the number of bosons in the condensate), 
we obtain: 



AS k = n^(k)[(P + l)i±^--^ 



(17) 



On the other hand, in the variational treatment of Eq. 
(Ilip , where we keep an accurate count of the number of 



bosons in the k = state as is done in the basis wavefunc- 
tions of Eq. (|12p . it can be shown 

.16] 

that the quantity 

Ai?k is given by: 



A£ k = n B v(k){ [C^(k 2 +3) + 1][- 



2c k 



, l-c k J 

(18) 

Using the expression of Ck given by Eq. ([SJ in Eq. p7|) 
above, one obtains the celebrated Bogoliubov spectrum 

0. Con- 



A-Ek = yk 2 (k 2 + 2), which is gapless as k 
versely, when the coefficients Ck of Eq. (|15j) are used in 
Eq. ([T8| . one obtains: 



A£ w 



n B v(k)^/Q 2 (Q 2 +2) 1 Q 2 ~ C^{k 2 + a), 



(19) 

which has a finite gap as k — > 0, Ai?k^o = 
gnsGj ^/cr((7 + 2Cd). For the values of n^a 3 , a and Cd 
considered in this paper, we obtain AEk-s-o = 0.98gns, 
which is comparable to the value of the gap gns predicted 
by the standard Hartree-Fock method, and by Girardeau 
and Arnowitt in Ref. [lii 

Conclusion. To summarize, in this paper, we have ar- 
gued that the decoupled way in which the Hamiltonian 
H = Ylu=io is diagonalized in the standard formula- 
tion of Bogoliubov's theory, where each and every mo- 
mentum contribution is diagonalized separately, is 
not appropriate. Diagonalizing the total Hamiltonian H 
directly leads to results that are markedly different from 
the results of Bogoliubov's method. More specifically, 
we find that the depletion of the condensate is smaller 
than what Bogoliubov's theory predicts, and that the en- 
ergy to excite a single boson from the condensate to the 
single-particle state with wavevector k has a finite gap 
as k — > 0. A more thorough analysis detailing further 
evidence in support of the above conclusions, including a 
more detailed discussion of the elementary excitations of 
the full Hamiltonian H, will be presented elsewhere [l6j |. 
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